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1. Introduction 

Recently, Hofava constructed an intriguing theory of gravity in four dimensions [1]. 
Motivated by models for condensed matter systems that exhibit anisotropic scaling 
phenomena, this Hofava-Lifshitz gravity (HLG) lacks four-dimensional Lorentz invari- 
ance but seemingly has better UV properties than Einstein gravity. A neat aspect of 
HLG is that for a particular choice of parameters the theory has a classical Weyl gauge 
symmetry while remaining second order in the time derivatives. This leads to a natural 
question: is the Weyl invariance anomalous? 
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In this note we study a technically simpler but related question: a matter system 
exhibiting appropriate anisotropic scaling may be coupled to HLG in a Weyl-invariant 
fashion; is there a Weyl anomaly for the matter system in the HLG background? We 
study the problem for a z = 3 free Lifshitz scalar in d — 4 and find an anomaly. 

To describe our result, it is useful to recall the work of [2] on the classification of 
Weyl anomalies in Lorentz-invariant theories in even dimensions. Consider the renor- 
malized effective action W[g; fi] obtained by integrating out the matter fields. Here g is 
the background metric and \i is a mass scale which may arise in removing logarithmic 
divergences. The Weyl anomaly is given by 



tp=0 



Following the terminology of [2], dilatation-invariant (i.e. /i-independent) terms in W 
give rise to "type A" terms in A, while //-dependent terms give rise to "type B" anomaly 
terms. Recall that in two dimensions there is only a type A anomaly, while in four 
dimensions the result is 

A\ d=4 = -a(Euler) + c(Weyl) 2 + bUR. (1.2) 

The first of these is of type A, the second is the type B anomaly, and the third may be 
removed by a local counter-term. 

These anomalies may be uncovered by studying n-point functions of the energy 
momentum tensor in flat space. The d = 2 case is a text-book example (see, e.g. 
[3]), where the study of the two-point function (T^a^T^O)) quickly leads to the 
anomaly. In four dimensions, since the non-trivial terms on the right-hand side of ( |1.2|) 
are quadratic in curvature, a study of three-point functions is necessary to directly 
determine a and c. However, the scale dependence of W[g; fj] is apparent at the level of 
the two-point functions. To see this, note that the scale dependent term has a schematic 
form 

W\g;n] =cjd 4 x v ^Clog(n//i 2 )C, (1.3) 
where C is the Weyl tensor and □ is the covariant Laplace operator. It follows that 



d_ 
dfi 



H— (T^(x)T pCT (0)) = cA^J^ix), (1.4) 



where A^ pa is a fourth-order derivative operator consistent with = 0, d^T^ = 
and Bose symmetry. Indeed, as shown in [4], a naive position-space computation of the 
two-point function yields a singular distribution, and the expected scale dependence 
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may be uncovered by carefully regularizing this distribution. Thus, the coefficient of 
the B-type anomaly may already be extracted from the two-point function. 

We have studied the two-point functions of the conserved quantities for the z = 3 
Lifshitz scalar. We find that while no type A terms arise at this order in the background, 
a logarithmic divergence leads to a scale dependence of the two-point functions and thus 
to a type B anomaly. 

Our result, while perhaps not very surprising, merits attention for several reasons. 
First, it bears on the UV behavior of HLG theories coupled to matter. As observed 
in [1], to achieve the desired UV improvements, it is necessary to work at points of 
enhanced gauge symmetry. Classically, the Weyl-invariant point is one such choice, 
and our results rule it out quantum mechanically once HLG is coupled to a scalar. 
Second, our work is an illustration of the complexity associated to even the simplest 
computations in anisotropic theories of gravity, and we hope that the techniques de- 
veloped herein may be of use in further investigations. Finally, the Weyl anomaly in 
a Lorentz-invariant theory is well-known to contain important physical information; 
our work may be viewed as a first step in learning how to quantify and extract this 
information in anisotropic theories. 

The rest of the paper is organized as follows. In section |2], we introduce the z = 3 
free Lifshitz scalar, compute its propagator and couple it to the Weyl-invariant HLG 
theory. In section [| we develop some position-space techniques and test these in the 
familiar relativistic d = 2 and d — 4 examples. Finally, in section ^, we apply these 
techniques to 2-point functions of conserved quantities in the z = 3 theory and find 
the scale dependence symptomatic of a type B anomaly. We conclude with a brief 
discussion of our results and some observations on possible type A anomaly terms that 
follow from Wess-Zumino consistency conditions. 

2. The Lifshitz scalar in a gravitational background 

2.1 The free scalar in flat space 

Our starting point is the free scalar with dynamical exponent z in d = 4. The Euclidean 
action is 





where d 2 = <9j<9j is the spatial Laplacian, and = d T (f>. The propagator for <j) has a 
simple Fourier space representation. Defining A( solution to 




(2.2) 
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we see that 



. , . f dud 



dud 3 k e iWT+ikx 



For our purposes, an explicit position-space representation will be useful. Performing 
the uj and angular integrations leads to 



A(t,x) 



47r 2 ; 



3-2 Z _ smivzu? ) 
dv v * e j — j — . 

V*U2 



where u = x 2 (t 2 ) x l z . Performing the remaining integral, we have 



A(t,x) 



4:% 2 z ^ T(2n + 2) 

n=0 v ' 



-u) n . 



(2.4) 



(2.5) 



Plugging in z — 1, we find the usual propagator for a relativistic massless particle. Our 
interest is in the more exotic limit of z — 3, which leads to 



A(t,x) 



1 



z=3 



127T 2 



G{u)- l -\og{r 2 m 2 ) 



(2.6) 



where m is a scale introduced to make sense of the logarithm and to absorb an infinite 
constant; u = x 2 (r 2 )~ 1 / 3 , and 



G(u) = 



r(f) 

^r(2n + 2) 



(-«)». 



(2.7) 



G(u) is an analytic function in the complex plane. As we will see, the relative complexity 
of A(r, x) at z = 3 makes even free field theory computations a bit more laborious than 
in the familiar z — 1 case. We also note that z = 3 is a natural limiting case, in which 
(f> has zero scaling dimension and A(r, x) acquires the log r 2 m 2 term. 

2.2 The scalar in a curved HLG background 

The z = 3 Lifshitz scalar has a natural coupling to z = 3 HLG; moreover, it is easy to 
construct a Weyl-invariant coupling. To describe this, we will first review some basic 
features of z = 3 HLG. We follow [1]. 

The HLG theory is defined on a four-dimensional space-time equipped with a co- 
dimension one foliation, where the latter structure encodes the privileged role of the 
time direction. The degrees of freedom are familiar from the ADM decomposition in 
Einstein gravity: there is a metric on spatial slices, gij(r,x), a shift one-form iVj(r, x), 
and a lapse function v(t,x). The action is given by 



/ 



S'hlg = drd x ^fgv 



2 K 2 

— {KijICi - \K 2 ) + — CwC^ 



K 



2w 4 



(2.8) 
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where k, w and A are undetermined parameters, Cij is the Cotton tensor constructed 
from gij, and 

K tj = g tj - ViN 3 - VjNi (2.9) 

is the extrinsic curvature. The connection V is the Levi-Civita connection associated 
to g. 

Sklg is invariant under spatial diffeomorphisms, time reparametrizations, and for 
A = 1/3 under local Weyl rescaling. Denoting the infinitesimal parameters for these 
transformations by £ l (r,x), f(r), and u>(t,x), the action on the fields is 

Sg = C^g + fd T g + 2wg, 
5N = C^N + d T (fN) + d T ^g + 2ooN, 

Su = C^ + d T (fu) +3cuu. (2.10) 

Here denotes the Lie derivative with respect to the vector field £, and l denotes 
the interior product, viz. (^l^)j = ^9ji- ^ is a simple matter to show that the algebra 
closes. 

Having described the gravitational theory we have in mind, we now turn back to 
the z = 3 scalar. An action in the curved background that reduces to the correct flat 
space limit and has diffeomorphism and time reparametrization invariances is given by 

S = J dT^Xv^l^-^W^-^V 2 ) 3 ^}. (2.11) 

S is invariant, provided transforms as a scalar: 

5<j> = ?di<j> + fj>. (2.12) 

In fact, the kinetic term is also Weyl invariant if <p has (as expected) Weyl weight zero. 
So, to construct a Weyl-invariant action we just need to make a suitable modification of 
the potential term. This is easily achieved by introducing a Weyl-covariant derivative 
V. 

To describe the construction, let V denote the Levi-Civita connection associated 
to the Weyl-invariant metric g = v~ 2 ^g. Given a tensor T with Weyl weight q, we 
define VT via 

VT = VT-gA®T, (2.13) 

where A = ^dlogu. It is easy to see that under a Weyl transformation = qujVT. 

Furthermore, since dA = 0, the curvature associated to V is just the Riemann curvature 
associated to the metric g. Finally, two useful properties of V are 

Vi/ = 0, Vg = 0. (2.14) 
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Using the connection V we can construct a Weyl-invariant potential term. In fact, 
many choices are possible. In what follows, we will restrict our considerations to a 
two-parameter family of terms, replacing z/0(V 2 ) 3 with 



Sv = t;J drd3x \f9 p [aViVj-V^VVvV + (3V 2 V k <pV 2 V k <p + 7 V fc V 2 0V fc V 2 < 

(2.15) 

with non- negative parameters a, f3, 7 constrained by a + f3 + j — 1. Note that we have 
chosen these terms to keep the improved action positive-definite. Such a requirement, 
while not very sensible for theories coupled to Einstein gravity, does seem to make sense 
in the context of HLG, since the latter minimal action is itself positive-definite. 



2.3 Conservation laws in flat space 

Having derived the curved-space action, we can follow the usual logic to find the con- 
servation laws in the flat background (i.e. N — 0, v — 1, = 5ij). To this end, we 
compute 



with 1 



5S\ 



flat 



drd x 



-bvE - -5g u T u - 5N^ 



P; 



(2.16) 
(2.17) 



e = ±0- + f r ijk + 



(2.18) 



and 



-a(j) ijk (j) ljk + (P + ^)d 2 (j)id 2 (t)i + (2a - ^)<\>u k d 2 <\) k + 2a<p iljk <f) jk - 7<fe(<9 2 ) 2 

-(a + ^d 2 <Pad 2 <P + Md 2 ) 2 <Pi + Md 2 ) 2 <f>i - a(<p ik d 2 <p lk + <\> lk d 2 <\> ik ) 

-(a + 7 )<9 2 iifc fc . (2.19) 



By using the variations in ( p.lOj ), we extract the following conservation laws. Time 
reparametrization invariance leads to 



E = / d 3 xd T £ = 0; 



(2.20) 



1 We will use a short-hand: 4>i 1 ---i n = di 1 ■ ■ ■ 
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spatial diffeomorphism invariance leads to 




(2.21) 



and, finally, Weyl invariance implies 



W = 3S + Tu = 0. 



(2.22) 



A short computation shows W = identically, while E — and J, 



m 



hold modulo 



the equation of motion. 

2.4 Quantum corrections and the anomaly 

Having described the classical conservation laws in flat space, we have two possible 
routes to determining the form of the anomaly. The first, which would produce the 
most satisfactory results, would be to construct the renormalized quantities p, 
and £ in an arbitrary gravitational background. This sort of computation is already 
quite heroic in relativistic theories [5] and would be challenging to implement here. 
The second, which is how the relativistic conformal anomaly was discovered in the first 
place [6], is to study the n-point functions of Pj, and £ in flat space and extract the 
presence of the anomaly from these. It is this second method we will pursue below. 

In order to preserve the conservation laws J m = 0, E = 0, and W = we must 
find local counter-terms so that, for example, 



In Lorentz-invariant theories this would be a simple one-loop computation, easily 
performed in either momentum or position space. In case of a Lifshitz theory, the 
problem is more involved. Some methodology is available for handling loop integrals 
in Lifshitz theories [7], but the expressions for even simple diagrams with non-zero 
external momenta are quite forbidding. The expressions simplify for one-loop graphs 
with zero external momenta, leading to, for instance, tractable computations of beta 
functions and gap equations [8-10]. To compute correlators at non-zero external mo- 
menta, we found it easier to consider the computation in position space. Setting up 
this computation in a convenient regularization scheme will occupy us in what follows. 

3. A regularization and position space computations 

In this section we present a regularization scheme and apply it to anomaly computations 
for relativistic scalars, where the results are well-known. 



((P + d J T ji )(x)£(0)) + d T (P(x)£(0)) ct + d^Tj&WO))* = 0. 



(2.23) 
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3.1 The regularization scheme 

In the course of position-space computations in a relativistic, rf-dimensional, scale- 
invariant theory, one typically encounters singular distributions such as |a;|~ 2d . A nice 
way to deal with these distributions is via the method of differential regularization [4, 
11]. However, since our final goal is to study the anisotropic theory, we will instead 
introduce a smeared propagator as was considered in [12]. That is, we introduce a scale 
e and replace the standard propagator A(x 2 ) with 

A e (x 2 ) = A(x 2 + e 2 ). (3.1) 

This leads to a representation of the Dirac 5-function via Green's equation, 

Se(x) = -d 2 A e (x 2 ) , (3.2) 

since 

lim / d d x S e {x)f(x) = /(0) (3.3) 

e->0j 

for any bounded function f(x). 

In our computations we will encounter distributions D e (x) that satisfy two basic 
properties: D e (x) is a smooth bounded function with lim e ^ D e (x) = for any i ^ 0; 
and D e (x) ~ e fc |rr|~™ for large \x\, with k > and n > d. We can find a convenient 
representation for such distributions by integrating them against a smooth, bounded 
test function f(x). Specifically, we have 



n—d—k 

e 



J d d x D e (x)f(x) = J d d y D 1 (y)e k - n+d f(ey) 
= J d d yD 1 (y)e k - n+d 

+ J d^DMt 



r k—n+d 



ml 

m=0 



(3.4) 



Since f(x) is smooth and bounded and J d d yDi(y)\y\ m < oo for m < n — d, it follows 
that the second line is a convergent integral for e ^ 0. Furthermore, by Taylor's 
theorem the integrand in the second line scales as e for small e, so that this remainder 
term vanishes in the e — > limit. Keeping the terms that do not vanish as e — > 0, we 
obtain a representation for the distribution: 

n-d-k m +k+d-n ( _\m 

- E r-^-S^-d^J^ix), (3.5) 

771. 

m=0 
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where the coefficients S 11 '" 1 " 1 are obtained by computing the convergent integrals 

S n- lm= f d d yD 1 (y)y^---y l ™. (3.6) 



This regularization may be adapted to the Lifshitz scalar with two noteworthy 
modifications. First, a look at the z = 3 propagator given in (|2.6|) shows that it is 
sufficient to do the smearing in the time direction. That is, we will replace 

A(r 2 , x 2 ) -> A e (r 2 , x 2 ) = A(r 2 + e 6 , x 2 ). (3.7) 

Second, while integrating the distributions we will encounter against a smooth bounded 
test function still leads to sensible representations, the integrals that must be evalu- 
ated are quite a bit more difficult, involving a large number of terms of products of 
hypergeometric functions. We have opted to handle these numerically. 

3.2 The relativistic scalar in two dimensions 

As a warm up exercise we treat the well-known case of the relativistic scalar in two 
dimensions using the prescription of e regularization introduced above. 
The action of a free scalar in two dimensions is 

S = J d 2 xd^ct) , (3.8) 

where /x, u, ■ ■ ■ = 1,2 and we work in a Euclidean signature so indices are raised and 
lowered using the metric 5^ u . The regularized scalar propagator is 

A £ (x) = --^log(:r 2 + e 2 ) , (3.9) 
leading to a representation for the Dirac 5-function 

«■ = ' (3 ' 10) 

The energy-momentum tensor is 

= d^d v (j> - -S^dx^cj) . (3.11) 

This tensor is conserved up to the equations of motion and is identically traceless. 

The first step in the analysis of the anomaly is to check whether the symmetry under 
diffeomorphisms is violated quantum mechanically in the presence of a gravitational 
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background. This is captured to first order in the metric perturbation by the two-point 
function 

Cvap = (d»T^(x)T a(3 (0)) . (3.12) 
Performing the Wick contractions and using the regularized propagator we have 



Cvap = d a d 2 A e {x)d Pv A e (x) + d p d 2 A e {x)d au A e (x) - 5 afS d x d 2 A e (x)d uX A e (x 



2e 2 



ir 2 (x 2 + e 2 ) 4 



x 2 — e 2 
x 2 + e 2 



5 a pX u — 



XfyX RX V 

xr + e z 



(3.13) 



Note that lim e ^ Cuaplx^o = 0, so as expected, the violation of the energy-momentum 
conservation is local. In order to extract the local contact-terms the procedure outlined 
in subsection |3.1| is used. One then has 

d 2 y 



d xC„ a p(x)f(x) 



7T 



\y 2 + iy 



y 2 -i 



y 2 + 1 



y 2 + l 



y a ypyu 



(3.14) 



where we have changed the integration variable to x = ey. Expanding f(ey) in powers 
of e and doing the angular integration (keeping in mind that terms with odd powers of 
y vanish in the integration due to the spherical symmetry) we finally get 



a 



vafi 



07re z 4o7T /47T 



vafJ 



X 



— ^{8fMxS u /3 + 8^b va — 8^ u S a p) 



(3.15) 



where only the terms which do not vanish in the limit e — > have been kept. Hence, 
we see that energy-momentum conservation is violated quantum mechanically. 

In order to preserve energy-momentum conservation (and the related symmetry 
under diffeomorphism) we introduce the following local counter-term 



6®(x) , 
(3.16) 



where On = tn V d v ■ Restoration of energy- momentum conservation requires that 



48vr ' 67re 2 
However, this clashes with the Weyl symmetry since 



Mx)T a/3 (0)> 



ct 



12tt 



' ~dj^ 2 \x) 



(3.17) 



(3.18) 
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Hence, the regularization scheme used here reproduces the familiar Weyl anomaly (e.g. 
in [3]). 

Finally, in order to determine the type of the anomaly we follow [4] and compute 
the dependence of the two-point function on the regularization scale e 

e^{T,„(x)T pa (0)) . (3.19) 
Treating the resulting distribution using the regularization scheme, we find 



e§- e (T»„(x)T p!7 (0)) = i- 



_ 2 p _ 

£ 2 r P v P° 



+ I pupa (-.(.Ippva Ipavp) 



5 (2 \x) + 0(e) , 

Ipvpa = fipi/dpa + dpo-dpi; ■ (3.20) 



At first sight, it seems as if the logarithmic derivative has a finite term and a scale 
fi needs to be introduced to properly define (T^ u T pa ). However, the indices in two 
dimensions can take only two values, so the finite term actually vanishes and 

e^-((T^(x)T p(T (0)) + (T, u (x)T pa (0)) ct ) = . (3.21) 
As expected, the anomaly is type A in two dimensions. 
3.3 The relativistic scalar in four dimensions 

In order to demonstrate the regularization method applied to a known case when there 
is a type B Weyl anomaly, we consider a scalar field coupled conformally to gravity. 
The regulated propagator in this case is 

A <« = STO ' (3 ' 22) 

and the regularized version of the 5-function satisfying d 2 A t (x) = —S t (x) is given by 



«■<*) = ' (3 ' 23) 



The improved energy-momentum tensor is [4] 

1 



Tp U = d^dA - + ■ (3-24) 
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The first thing to be done is to compute the violation of the energy-momentum 
conservation in the two-point function. Using the same techniques as in the d = 2 
example, we find 



(d»T^(x)T pa (0)) = 



1 



60tt 2 



1 



AiSpvdv + A 2 (5 up d tT + 5 ua d p ) + 



+ ^7 



Mdvpa + A±b pcj d v d 2 + A 5 (S up d a + 5 ua d p )d 2 



+ 



+ [A 6 d pai/ + A 7 5 pa d u d 2 + A 8 (5 up d a + 5 va d p )d 2 ] <9 2 }^ 4 )(x) , (3.25) 



where the coefficients are 



A 1 = A 2 = 1 , A 3 



1 1 

"6' 4 " 12' 



. 1 . 7 1 

^""24' ^-"192' A8 "64- 



(3.26) 



Next, we turn to the violation of the Weyl symmetry in the two-point function. 
Here the four-dimensional theory differs from the two-dimensional one since the trace 
only vanishes up to the equation of motion. The most general possible contact term 
allowed by Lorentz and Bose symmetries and dimensional arguments is 



(TH(x)T pa (0)) 



60tt 2 



\B x b pa + \{B 2 5 pa d 2 + B 3 d p d a ) + (B 4 5 pa d 2 + 



+ B b d p d a )d 2 



5 (4) ( 



X) 



(3.27) 



Direct computation leads to the coefficients 



B l = 6 , B 2 = 



1 



B, = — 



B, 



B, = — 



1 



(3.28) 



11 

12 ' " J 6 ' ~96 ' " J 24 

On the other hand, the most general local counter-term consistent with Lorentz 
and Bose symmetries and dimensional analysis is 

1 f 1 

(T pu (x)T pa (0)) ct = —± [Ci W + CziSupdva + S pa 5 up )] + 

+ CeiS^ + 5 M( A P )<9 2 ] + CVVM^ 2 ) 2 + 

• <"sl'V/''V" ' ) + C^Ip Vpa d I ('\[){l lipl •■prri'p 

+ C 11 dp l/ p a \s^(x) . 



)d 2 + 
(3.29) 
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= c 2 = 


-1 , 


c 3 = - 


1 

"12 ' 


c 6 


1 

~ ~4 ' 


c 7 = 


11 
192 


C u 
2 


c 9 


1 

~48 


2 


, C10 


1 

~ 32 



By requiring these counter-terms to cancel the contact terms in (|3.25|) and (|3.27 ), we 
get a one-parameter family of solutions: 

Ca — — , Ck = — , 
8 6 

3 3Cn 

■ ^ • (3.30) 

So after adding the counter-terms the energy-momentum tensor is conserved and trace- 
less in the two-point function. This matches the well-known result that the trace is 
non-vanishing only in the three-point function. 

However, by considering the logarithmic derivative 

e^-(T^(x)T pa {0)) = g^j^j [DiS^ + D 2 {5„ p 5„ a + S^S vp )] + ^ [D 3 5^5 pa d 2 + 

+ D^S^pdva + 5 pa 5 up )d 2 + D^I pvp(T + D§{I ppva + I vppa )\ + 
+ D 7 5 p J pa (d 2 ) 2 + DstfwS** + <W<W)(d 2 ) 2 + D 9 I^ pa d 2 + 

+ Ao(W + W)^ 2 + D n d^8®{x) (3.31) 

it is possible to determine whether a scale signaling the violation of the Weyl symmetry 
is introduced. Direct computation of the coefficients using the regularization and the 
techniques presented above yields 



Di = D 2 = -4 , D 3 = \, Di = ~, D 5 = ~, As = 7 , 

3 2 6 4 

I „ I „ 1 „ 1 1 



D 7 = — , D 8 = , D 9 = , D l0 = — , Du = . (3.32) 

7 24 ' 8 16 24 ' 16 12 v ; 

The logarithmic derivative of the renormalized energy-momentum two-point function 
(which includes the contribution for the counter-terms) is 

e^{(T, v {x)T P M) + (T M ,(x)T pCT (0)) ct ) = —L-^[25, u 5 p(7 (d 2 ) 2 ~ 

lupoid ) 9 

-4^ p(T ](5 (4) (x) . (3.33) 

We see that the terms which diverge as e — > are canceled by the counter-terms and 
only the finite piece remains. This result matches the one in [4] up to a minus sign 
coming from the derivative being with respect to a length scale and not a mass scale - 
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lending credence to our regularization method. Integrating the above equation leads 
to a logarithmic term log(/ie) with fi being a scale introduced in order to make the 
argument of the logarithm dimensionless. Thus the renormalized two-point function 
and hence the effective action contain a scale /i. This signals that the Weyl symmetry 
is broken by a type B anomaly. 



4. The Weyl anomaly 

In this section we show that the conformally-coupled scalar in a Hofava-Lifshitz back- 
ground induces a Weyl anomaly. 

The propagator of a scalar coupled conformally to a flat Hofava-Lifshitz background 
( |3.7| ), together with the Green's equation DA e (r, x) = —S e (r,x), where □ = c% + (<9 2 ) 3 , 
leads to the regularized ^-function 

^ X) = 18^ + e 6 ) 2 + 3 ) F(Ue) ' 

oo V (2n ,i\ 

n») S £ 1 ^ TT yf HO", (4.D 

n=0 V ' 

where u e = ^ r2 ^6y/3 ■ We demonstrate in Appendix |A] that 

/**,«r,,)-l. (4.2) 

4.1 Counter-term analysis 

As in the relativistic examples, we begin with the quantum violation of energy and 
momentum conservation Ward identities. As expected from dimensional analysis, ro- 
tational invariance, parity and time reversal symmetry, they are of the form 

(E(t)£(0, 0)> = {^Assi + A g£2 d^j d T 5(r) , (4.3) 
(E(t)TI(O, 0)> = 5 ij (^Aen + A £T2 dA d T 6(r) , (4.4) 



{J\r,x)Pi (0,0)) 



(£(r)P*(0,0)> = , (4.5) 

1 



2 A JP1 5 ij d T + (A JP2 d» + A JP ^d 2 )d T 



S^(t,x) , (4.6) 
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(.r(T,x)T lm (0,0)) = { (S mi S kt + S ll S mk )d k 



-A JT1 + -^A JT2 «9 2 + \A JT3 {d 2 ) 2 + 



+ A JTi (dy + A JT5 d: 



+ 5 lm d l 



AjT6 + 



+ - 4 A JT7 d 2 + - 2 A JT% {d 2 ) 2 + A JT9 (d 2 ) 3 
e 4 e 2 

+ A JTW d 2 



+ A JT13 (d' 



2\2 



_|_ Qilm 



e 4 



(4.7) 



(J m (r,x)£(0,0)) = d r ' 



- 6 A J£1 + - 4 A J£2 d 2 + - 2 A J£3 (d 2 ) 2 + A J£i (d 2 ) 3 + 
e° e 4 e z 



+ A J£5 d 2 



The Weyl Ward identity is not corrected before the introduction of counter-terms since 
it vanishes identically. 

We now illustrate the procedure for computing the coefficients in the contact terms 
by evaluating Ajpi, Ajp 2 and Ajp 3 in (|4.(j|). The Wick contractions yield 



C ij = (J 1 (t, x)P j (0, 0)) = -d*d T A e (r, x)d j 5 e {r, x) - d lj A e (r, x)d T 5 e (r, x) . (4.9) 
Substituting the expressions for A e and 5 t , it can be verified that this distribution meets 



the requirements listed in subsection ^Tj and the use of the test function approach is 
possible. 

To extract the expression for the distribution we consider 



drd 3 xC ri (T,x)f(r,x) 



(4.10) 



for a smooth and bounded test function f(r,x). We first change variables to x % = eu l , 
r = e 3 t, and expand the test function f(eu, e 3 t) in powers of e. Performing the angular 
integration and changing variables to v — u 2 /(t 2 + 1) 1//3 in order to disentangle the u 
and t integrations, we obtain / = I\ + I 2 , where 



h 
h 



277r 3 e : 



Af8 r 



dt t' 



547T 3 



d T dkif 



,o (t 2 + 1) 

dtt 2 
(t 2 + l) 5 / 2 



— & J™ dun 1 ' 2 (^-uH^u) + H 2 (ufj 



—(5 ri 5 kl + 25 rk 5 u )J 1 + -5 ri 5 kl J 2 
15 3 



(4.11) 



- 15 - 



H^u) = (G' + 2uG")F[ + 6G"F 1 , 
tf 2 ( u ) = uG ! F[ + QG'Fi , 



(4.12) 



and 



J 2 



duu 5/2 H!\ 



u 



duu 3/2 H 2 (u) . 



(4.13) 



Ji and J 2 have been evaluated numerically and are convincingly given by 

T 1 T T 27 



Hence, 



A JP1 = 2.02 x 10" 



.4 



JP2 



1 

'720tt 2 ' 



.4 



JP3 



(4.14) 



(4.15) 



1607T 2 

The rest of the coefficients are determined by the same procedure. The numeric results 
are given in Appendix |B|. 

The most general local counter-terms with the appropriate dimension consistent 
with three-dimensional rotational symmetries, parity, time reversal and Bose symmetry 
are 



(££) 



ct 



-Cfici + - d C ££2 d 2 + \c ££3 (d 2 ) 2 + C ££4 (d 2 ) 3 + 



5 (4) (r,x) 



(4.16) 



(ST ij ) 



ct 



5 y 



-zCgTl + ~jC £ T2d + -^C £ T?,(d ) + C £ T4(9 ) + C £ T<od r 



—tC £ t& + -^C £ T-jd 2 + C £ Ts(d 2 ) 2 



(£P%t = C £P1 d i d T 5^(r,x) 



6®( 



T,X) 



ct 



^ ( ^C PP1 + C PP2 9 2 ] + Cpp 3 ^ 



5 (4) ( 



[P l T ]k ) ct = [6> k C PT1 d i a r + (S ij d k + 5 ik d ] )C PT2 d T ] 5 (4) (r, 



+ 

(4.17) 

(4.18) 
(4.19) 
(4.20) 
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/rpijrpkh 



Ct 



5 ij 5 kl 

ikx.il i xilxik 



+ 



+ (W + 8 U S J 
+ Ctt9(^ 2 ) 3 + Cttio^t 



-^Ctt& + -^CttiQ 2 + -^CTTs{d 2 ) 2 + 
e° e 4 e z 



-Cttii + -^Ctti29 2 + Ctt13(9 



,2\2 



+ 



-Cttw + 



2\2 



^ (4) (r, 



(4.21) 



The counter-terms can restore both diffeomorphism and Weyl symmetries if the coef- 
ficients of the quantum corrections satisfy 





3^4 j£5 + 3A JT10 + 2A JT5 


= o , 


3A J£4 


+ ^4jti3 + 2Aj T4 + 3A JT9 


= o , 




Aes2 + ^4fT2 


= o , 




3A J£1 + 2A JT1 + 3A JT6 


= o , 




Asti — Ajsi 


= o , 




3A ££X — 2A JT i — 3Aj T6 


= o , 


3A J£2 


+ ^4jtii + 2A JT2 + 3Ajt7 


= o , 


3A J£3 


+ A JTU + 2Aj T3 + 3Aj T8 


= . 



(4.22) 



The numerically-computed values of the coefficients satisfy the relations ( |4.22| ) - 
the Ward identities are preserved at the two-point function level. 

4.2 Scale dependence 

We turn our attention to the derivatives of the two-point functions with respect to 
the regularization scale. These are local, so it makes sense to use the test function 
procedure. 

The derivatives of the correlation functions 



oe 



,)'■' ( + A^tf ) ■■■■■A% A 0'0> 



5^(x,t), (4.23) 



5 (4) (x,r) 



(4.24) 
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are indeed non-zero and in particular Ap g PZ is independent of a, (3 and 7 (see Appendix 
H). Therefore, the correlation functions must have a log(/xe) term regardless of the 
values of the couplings, signaling the existence of a type B Weyl anomaly in the theory. 



5. Discussion 

We have shown through an explicit computation that the z = 3 Lifshitz scalar has a 
type B Weyl anomaly when coupled to a background of Hofava-Lifshitz gravity intro- 
duced in [1]. We suspect a similar result will hold in pure HLG itself. Experience with 
relativistic theories would suggest that exact Weyl invariance only comes at the price 
of perturbative unitarity. 2 It is conceivable, though perhaps unlikely, that anisotropic 
gravity might evade the positivity requirements. In that case, a modification of HLG in- 
volving additional degrees of freedom and gauge symmetries (such as might follow from 
considering some anisotropic locally supersymmetric theory) could be Weyl invariant. 
Such developments might be interesting, but the idea is uncomfortably reminiscent of 
the familiar "Stone Soup" tale. 

The computational technique we use is difficult to extend to three-point functions, 
which would be required to find a type A anomaly. It would be useful to develop 
more powerful techniques for flat space computations. A more ambitious and difficult 
undertaking would be to generalize the point-splitting techniques reviewed in [5] to 
the anisotropic case and thereby compute the renormalized Tij, S and Pi in an arbi- 
trary HLG background. Together with results on chiral anomalies, such as obtained 
in [8], these will yield important constraints on the structure of correlators of conserved 
currents in Lifshitz theories, analogous to those obtained in [4] for conformal theories. 

A natural step in studying the structure of anomalies is to consider the Wess- 
Zumino consistency conditions. 3 We will end our work with a brief look at the consis- 
tency condition for Weyl invariance. Consider the one-loop effective action W[g, N, u] 
and assume there exists a regulator that preserves diffeomorphism and time reparamet- 
rization invariances. On general grounds the Weyl variation of W is given by a local 
functional 



<LW = J drd 3 x y/guAcu, (5.1) 



2 It has been shown that conformal supergravity coupled to a certain super Yang-Mills theory has 
no Weyl anomaly [13]; however, that theory is haunted by the ghosts of conformal graivty. Moreover, 



it has been argued in [14] that unitarity implies positivity of the anomaly coefficient a in (5.3) 
3 These are concisely reviewed in the context of Weyl anomalies in [15]. 
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where the local function A transforms as a scalar under diffeomorphisms and time- 
reparametrizations and satisfies the Wess-Zumino consistency condition: 

[<L 15 <L 2 ]W = J drd 3 x ^fgv [u; 2 <Li - ^i<L 2 ] -A = 0. (5.2) 

In Lorentz-invariant theories this is a constraining requirement. For example, in d = 4 
it allows just three independent purely gravitational parity-invariant terms in A [16]: 

A = -a(Euler) + c(Weyl) 2 + bUR. (5.3) 

In the anisotropic HLG theory many more terms are possible. Let Rij be the Ricci 
tensor constructed from the Weyl-invariant metric g = v~ 2 / 3 g. Clearly, any scalar A of 
the schematic form 

A = VRVR + RVVR + RRR, (5.4) 

where the indices are contracted with the metric g, will trivially satisfy the consistency 
condition. These include the square of the Cotton tensor but clearly contain additional 
terms. It would be interesting to classify these and to determine which, if any, may be 
eliminated by local counter-terms. 
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A. The regularized anisotropic ^-function 

We show here that (|Q|) satisfies 

I = J drd 3 xS e (T,x) = 1, (A.l) 

as required by a representation of Dirac's 5-function. Changing variables x = eu, t = e 3 t 
and doing the angular integration and then changing variables again v = u 2 (t 2 + l) -1 / 3 
to disentangle the t and v variables yields 



I = — 

9?r 



1 f°° dt f°° , r- ( \ . . . A , 
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Plugging in the result of the t integration and integrating by parts in v, we have 

I = — / dvv 1/2 F(v) = — / dvF(v 2/3 ) . (A.3) 
3vr Jo 9tt y 

Fortunately, F(t> 2//3 ) has a nice integral representation 



, /"°° sin [y^Vl 

F (v 2/3 ) = lim / dy ^5 ^ . A.4 

p-+WJo y 1/3 w p 

Interchanging the v and y integrations and using the integral 

00 , sin(axf) v / ^2( 1 - 2 P)/fa- 1 /Pr(^) 

cte — i - = ^- , p > 1 , a > A.5 

to analytically continue to p = 1/3, one finally finds 7=1. 
B. Counter-term coefficients 

In this section we give the results of numeric computations for the contact term coeffi- 
cients arising in section ^. In the coefficients that follow we have extracted an over-all 
factor of 10~ 3 . 

^.££1 = 1-41, Af£2 = \Ag£i, Asti = —Assi, A£T2 = —jAssi] (B.l) 

A JP1 = 2.02, A JP2 = -0.14, A JP3 = 0.63; (B.2) 
Ajei = —A ££1 , A J£5 = ~A ££ i, 



(B.3) 



-0.28, 







A J£2 = -0.02 - 0.79a 


- 2.46/3 


- 1.357, 






A J£3 = 0.35 + 0.79a - 


1.40/3 - 


0.097, 






A J£4 = 0.30 + 1.24a - 


0.62/3 + 0.267; 


AjTl 


= 2 


.53, A JT2 = 0.22, 


AjT3 


= -0.04, 


AjTA 


= 


.27, A JT5 = 0.63, 


AjT6 - 


= —Ajtio 


AjT7 


= 


.02 - 0.70a + 1.81/3 + 


•147, 




AjT8 




■0.35 - 1.96a + 1.33/3 - 


0.64 7 , 




AjT9 




■0.30- 1.98a + 0.81/3 - 


0.517, 




Ajtu 


= 4. 


,03a + 1.52/3 + 3.197, 






AjTYl 


= 3 


,60a + 0.30/3 + 2.287, 






AjTIZ 


= 1 


,70a- 1.10/3 + 0.227. 







(B.4) 
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In the 



logarithmic derivatives we get 



A%, = -4.04 , 

App2 = 3.80 , 
/ii°g — _i 97 

A 1 ^ = -6.33 + 4.22(a + /3 + 7) - 6.33a 2 - 12.67a/3 - 6.33/? 2 - 80.21^7 + 

+ 4.22/37 - 56.997 2 , 
A%$ 2 = -4.47 - 2.83a + 3.07/3 - O.867 + 11.32a 2 + 8.89a/3 - 2.44/3 2 - 

- 78.27a7 + 28.38/37 - 65.5I7 2 , 

Aj?| 3 = -2.17 - 7.15a + 3.62/3 - 1.777 + 2.59a 2 + 8.33a/3 - 3.69/3 2 - 

- 54.66a7 + 14.58/37 - 43.82 7 2 , 

A£f 4 = -1.20 - 5.42a + 2.45/3 - 0.277 - 13.40a 2 + 2.37a/3 - 1.13/3 2 - 

- 31.02a7 + 3.54/37 - I6.8I7 2 , 

4°! 5 = • (B.5) 
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